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Abstract: We review recent results on the 9 dependence of the ground-state energy 
and spectrum of four-dimensional SU(iV) gauge theories, where 9 is the coefficient of 
the CP-violating topological term FF in the Lagrangian. In particular, we discuss 
the results obtained by Monte Carlo simulations of the lattice formulation of QCD, 
which allow the investigation of 9 dependence around 9 = by determining the 
moments of the topological charge distribution, and their correlations with other 
observables. The results for N = 3 and larger values of N support the scenario 
obtained by general large- N scaling arguments. 
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1. 6 dependence of the ground-state energy 

Four- dimensional SU(iV) gauge theories have a nontrivial dependence on the angle 
9 that appears in the Euclidean Lagrangian as 

Ce = \f; v {x)F; v {x) - l 9^- 2 e, vpu F^{x)F; rj {x) (1.1) 

2 

where q(x) = ■^2^fiupaF^(x)F pr7 (x) is the topological charge density. Semiclassical 
instanton solutions are one example of field configurations which have nontrivial 
topological properties. Moreover, the most plausible explanation of how the solution 
of the so-called U(l)^ problem can be compatible with the 1/N expansion (performed 
keeping g 2 N fixed [|TJ) requires a nontrivial 9 dependence of the ground-state energy 
density F(9) §,!, 

exp[-VF(0)] = J [dA}exp f- J d A xC^j (1.2) 

where V is the volume. Evidence for such a dependence has been obtained by exploit- 
ing the lattice formulation of the theory, using numerical Monte Carlo simulations, 
as will be described in Section ||. The complex nature of the 9 term in the Euclidean 
QCD Lagrangian makes the Monte Carlo studies of the 9 dependence quite hard, 
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since the lattice action corresponding to the Lagrangian (1.1) cannot be directly 
simulated for 9^0. Nevertheless, important information on the 9 dependence of 
relevant physical quantities, such as the ground-state energy and the spectrum, can 
also be inferred from results at 9 = 0, by expanding them about 9 = and computing 
the coefficients of the expansion [|], [|. The 9 dependence is particularly interesting 
in the large-iV limit where the issue may also be addressed by other approaches, 
such as AdS / CFT correspondence applied to nonsupersymmetric and non conformal 
theories, see e.g. Ref. ||. 

We introduce a scaling energy density 

fW = (1-3) 

where AF(9) = F(9) —F(0) and a is the string tension at 9 = 0. By expanding f(9) 
around 9 = 0, one can study its 9 dependence in the region of small 9 values. The 
function f(9) is conveniently parametrized as 

/(*) = 1 -C9 2 s(9), (1.4) 
where C is the ratio x/ " 2 an d X is t ne topological susceptibility at 9 = 0, 

X = J d*x(q(x)q(0)) = ^ (1.5) 

where Q = f d A xq(x). s{9) is a dimensionless function of 9 such that s(0) = 1. 
The function s{9) can be expanded around 9 = as 

8 (9) = 1 + b 2 9 2 + M 4 + • • • • (1.6) 

The coefficients of the expansion of f(9) are related to the zero- momentum n-point 
connected correlation functions of the topological charge density, and therefore to 
the moments of the probability distribution P(Q) of the topological charge Q. If 
s{9) = 1, and therefore hm = 0, the corresponding distribution P(Q) is Gaussian, 
i.e. 

i f Q 2 



p( ^wrH-#>J- (l7) 

Therefore the coefficients &2n of the expansion of s(9) parametrize the deviations from 
a simple Gaussian behavior. For example, he first non-trivial correction is given by 

k = -^, (1-8) 

X 4 = ^[(gVo-3((gVo) 2 ]. (i.9) 

It has been recently shown (see also 0) that correlation functions involving 
multiple zero-momentum insertions of the topological charge density can be defined 
in a nonambiguous, regularization-independent way, and therefore the expansion 
coefficients bm are well defined renormalization-group invariant quantities. 
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2. Behavior in the large- N limit 



Witten argued |§ that in the large- N limit F(9) is a multibranched function of the 
type 

'6 + 2<Kk s 



F(9) = iV 2 min fc H ( - j (2.1) 

which is periodic in 9, but not smooth since at some value of 9 there is a jump 
between two different branches. This issue was also discussed in Ref. |l(|. More 
recently, the conjecture was refined leading to a rather simple expression for 
AF(9) in the large- N limit, that is 

AF(9) = A mm k (9 + 2nk) 2 + 0(1 /N) . (2.2) 

In particular, for sufficiently small values of 9, i.e. \9\ < n, 

AF(9)=A9 2 + 0(1/N). (2.3) 

Thus possible 0(9 4 ) terms are expected to be depressed by powers of 1/N. This 
conjecture has been supported using arguments based on duality between large-iV 
gauge theories and string theory [[II]]. It has also been discussed in a field-theoretical 
framework in Ref. \T2 . 

The large- N behavior of the coefficients b 2n of the expansion of f(9) around 
9 = can be inferred by using general large- N scaling arguments applied to the 
Lagrangian ( |1 . 1|) . They indicate the ratio 9 = 9/N as the relevant quantity in the 
large- iV limit of the ground-state energy, and more generally of the spectrum of the 
theory. Then we expect 

f(9)=N 2 f(9 = 9/N), (2.4) 
f(9) = 1 -CJ 2 {1 + b 2 9 2 + b 4 9 4 + •••), (2.5) 

where is the large- N limit of the ratio C = xl° 2 ■ Comparing with Eq. (|1.4p , one 
derives 

C = C 00 + c 2 /iV 2 + ..., b 2i = b 2i /N 21 + (2.6) 

We recall that a nonzero value of is essential to provide an explanation to the 
11(1)^4 problem in the 't Hooft large-N limit, and can be related to the 7/ mass 0. [| 
through the relation 

Xoo = J -Q^ + 0(l/N). (2.7) 

The quantity b 2 also lends itself to a physical interpretation, being related to the 
7/ — 7/ elastic scattering amplitude B. 
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3. Results for the first few terms of the expansion around 
6 = of the ground-state energy 



The 9 dependence of SU(iV) gauge theories has been investigated by Monte Carlo 
simulations of their Wilson lattice formulation. The lattice action corresponding to 
the Lagrangian cannot be directly simulated for 9 ^ 0, by virtue of the complex 
nature of the 9 term. On the other hand, the coefficients h 2n in the expansion of 
the ground-state energy F(9) around 9 = can be accessed by determining the 
moments of the topological charge distribution at 9 = 0. They are dimensionless 
renormalizat ion-group invariant quantities, which should approach a constant in the 
continuum limit, with 0(a 2 ) scaling corrections (a is the lattice spacing). 

Computing quantities related to topology using lattice simulation techniques is 
not a simple task. In the case N = 3 several methods have been employed to de- 



0,1 



Cooling, 



termine the topological susceptibility, see e.g. Refs. 
geometrical, heating techniques have been used to address the problems caused by 
power-divergent additive contributions and multiplicative renormalizat ions in defini- 
tions of the topological susceptibility based on discretized versions of the topological 
charge density operator q(x). These methods have their drawbacks, since their sys- 
tematic errors are not under robust theoretical control. 

A substantial progress has been achieved after the introduction of the Neuberger 
overlap formulation [|39], |4(| of fermions, which represented a breakthrough for the 
lattice formulation of QCD. Overlap lattice fermions satisfy the Ginsparg- Wilson 
relation and therefore preserve an exact chiral symmetry | 4"2"|| . As a by product, 
the index of the overlap Dirac operator [ J3| provides a well-defined estimator for the 
topological charge |4(| [|4|, §, [7j , which can also be used in pure gauge theories. This 
method circumvents completely the problem of renormalization arising in bosonic 
approaches, even though at a much higher computational cost. Using these methods, 
the topological susceptibility of the pure SU(3) gauge theory has been investigated 
in Refs. [0-[^3j, finally obtaining the accurate estimate |k| x r o — 0.059(3) (r is 
the length scale defined in |54]|). This value corresponds to C = xjo 2 = 0.029(2) 
(using H o-VVo = 1.193(10)). 

It is important to note that the results obtained by the (less computer-power 
demanding) bosonic methods are substantially consistent, see e.g. Refs. K7], W2 



|3"B| , (|, [52], showing their effectiveness although they are supported by a weaker 
theoretical ground. For example, we mention the results: [27] C = 0.027(4), obtained 
using the heating method, ||] C = 0.0282(12), obtained using cooling, and the more 
recent result @ C = 0.0259(11) . 

For larger values of N, results have been obtained only by the cooling method so 
far |33], |, |36[], up to iV = 8. They fit well the expected large- N behavior: C = + 
c 2 /N 2 , providing an estimate of C^, and therefore of the topological susceptibility 
in the large-iV limit: = 0.0200(43) P§, C M = 0.0221(14) §, = 0.0248(18) 
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[j36|j (the latter was obtained using N < 8 and keeping a fixed). These results are 
in substantial agreement with the large- N relation ([2.7]) . We stress that the good 
agrement for N = 3 of the cooling method with the more rigorous overlap result make 
us quite confident on the reliability of results for higher values of N, since there are 
no arguments to suggest that this agreement could be spoiled with increasing N 
(actually there are reasons in favor of improved agreement JIS], |5B|). An independent 
determination of using other methods would be welcome. 

Higher moments of the topological charge distribution provide estimates of the 
coefficients b 2n of the expansion of the scaling energy density f(9), cf. Eqs. (|1.4|) and 
( |1.6|) . In particular b 2 can be estimated using formulae (|1.8| , [O^). There are a number 
of results at N = 3, obtained by different approaches: Ref. H used the cooling 



method, Ref. |57| used the heating technique to estimate additive and multiplicative 
renormalizations in zero-momentum correlations of lattice discretizations of q(x), 



and finally Ref. |58[ used the most rigorous and CPU intensive overlap method. The 
results reported in Table [l] are in good agreement, suggesting that the systematic 
errors of the various methods are sufficiently small. We mention that the fourth 
moment of the topological charge distribution has been numerically investigated also 



in Ref. [38[, without arriving at any definite conclusion. 



N Ref. method 



The results of Table [I] provide robust evidence that b 2 is nonzero, and therefore 
that there are deviations from a Gaussian distribution of the topological charge. 1 
However, b 2 turns out to be quite small, indeed \b 2 \ <C 1. Thus deviations from a 
simple Gaussian behavior are already small at N = 3. 



There are also estimates for larger values 
of N, see Table [I], but only using the cooling 
method. Again, given the agreement found 
at N = 3, higher N results should be suf- 
ficiently reliable. They appear to decrease 
consistently with the expectation from the 
large-iV scaling arguments, i.e. b 2 ~ b 2 /N 2 
with b 2 « -0.2. 



3 [ I J cooling 
|57| heating 
p3 overlap 



4 
6 



cooling 
cooling 



-0.023(7) 
-0.024(6) 
-0.025(9) 
-0.013(7) 
-0.01(2) 



Table 1: Results for the coefficient b 2 
of the expansion (lL 



We also mention that the analytical prop- 
erties at 6 = have been recently discussed 



and numerically checked in Ref. [37 



1 An apparently contradictory result has been reported in Refs. ^Tl [33| for the expected large- 
tip 

volume probability distribution P(Q), i.e. P{Q) = (27r(Q 2 ))- 1 / 2 e ^Q 1 ) [1 + 0{1/V)}. A purely 
gaussian behaviour would imply an exact quadratic form for /(#), and in particular a vanishing 62, 
thereby contradicting the assumption of a generic expansion of f{0). 
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Overall, these results support the scenario obtained by general large- N scaling 
arguments, which indicate 9 = 9/N as the relevant Lagrangian parameter in the 
large-iV expansion. They also show that N = 3 is already in the regime of the 
large- N behavior. For N > 3 the simple Gaussian form 

f{9) » l - x e 2 (3.1) 

provides a good approximation of the dependence on 9 for a relatively large range of 
values of 9 around 9 = 0. 



4. 6 dependence at finite temperature 

Another interesting issue concerns the behavior of topological properties at finite 
temperature, and in particular their change at the finite-temperature deconfining 
transition, which is first order for N > 3, see e.g. Ref. ^TJ and references therein. 



This issue has been investigated in a number of numerical works, see e.g. Refs. [21, 27 



30| , |62| , |60| , pq] , using different methods. They show that the topological properties, 
and in particular the topological susceptibility x> var Y ver Y little up to T < T c . 
They change across the transition, where \ shows a significant decrease. Then, at 
high temperature T 3> T c , where the instanton calculus is reliable, a rather different 
scenario emerges 



33 



Concerning the large- N behavior (investigated by performing simulations at var- 



ious values of N > 3 |60], |3q[), the results indicate that x has a nonvanishing large- N 
limit for T < T c , as at T = 0, and that the topological properties, and therefore 
F(9), remain substantially unchanged in the low-temperature phase, up to T c . On 
the other hand, above the deconfinement phase transition, for T > T c , \ shows a 
large suppression, hinting at a vanishing large- iV limit for T > T c . These results sup- 
port the hypothesis put forward in Ref. [63]: At large iV the topological properties 
in the high-temperature phase, for T > T c , are essentially determined by instantons 
that are exponentially suppressed, i.e. behave as e~ N , and therefore the topological 
susceptibility gets rapidly suppressed in the large- N limit. 



5. 6 dependence of the spectrum 

Another interesting issue concerns the 9 dependence of the spectrum of the theory. 
The analysis of the 9 dependence of the glueball spectrum using AdS/CFT suggests 
that the only effect of the 9 term in the leading large- iV limit on the lowest spin- zero 
glueball state is that this state becomes a mixed state of ++ and h glueballs, as a 
consequence of the fact that the 9 term breaks parity, but its mass does not change 
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Ref. Hfl presented an exploratory numerical study of the 9 dependence in the 
spectrum of SU(iV) gauge theories. Again numerical simulations of the Wilson lattice 
formulation were employed to investigate the 9 dependence of the string tension a (9) 
and the lowest glueball mass M(9). Around 9 = one can write 

a{6) = a (l + s 2 e 2 + ...), (5.1) 
M(9) = M(l+g 2 9 2 + ...) (5.2) 

where a and M are respectively the string tension and the ++ glueball mass at 
9 = 0. Then the coefficients of these expansions can be computed from appropriate 
correlators at 9 = 0. In particular, s 2 can be determined U from the large-t behavior 
of connected correlation functions of two Polyakov lines at distance t and the square 
topological charge, such as 

(Ap(t)Q 2 } e=0 - (Ap(t)} e=0 (Q 2 }e= (5.3) 

where 

A P (t)= ^TrPt( ;0)TrP( 2;i)a ; 2 ;t) ) (5.4) 

X\,X2 

P(xi,x 2 ',t) is the Polyakov line along the x% direction of size L, and Q is the topo- 
logical charge. Analogously, the 0(9 2 ) term of the glueball mass can be obtained 
from appropriate connected correlation functions of plaquette operators and Q 2 . The 
0{9 2 ) coefficients s 2 and g 2 are dimensionless scaling quantities, which should ap- 
proach a constant in the continuum limit, with 0(a 2 ) scaling corrections. 

Ref. || obtained the first estimates of s 2 and g 2 using the cooling method to 
determine the topological charge, and for N = 3, 4, 6 to also check their large- N 
behavior. The 0(9 2 ) terms in the expansion around 9 = of the spectrum of SU(iV) 
gauge theories are small for all N > 3, especially when dimensionless ratios are 
considered, such as Mj s/a and, for N > 3, the ratios of independent k strings. For 
example we mention the estimates s 2 = —0.08(1) and g 2 = —0.06(2) for N = 3. One 
may also consider the 9 dependence of the scaling ratio 

M (V) M n , Q 2 



a(9) V° 



l + c 2 9 2 + ...), (5.5) 

where c 2 = g 2 — s 2 /2, thus c 2 = —0.02(2) for N = 3. Moreover, the 0(9 2 ) corrections 
appear to decrease with increasing N, and the coefficients do not show evidence 
of convergence to a nonzero value. This is suggestive of a scenario in which the 9 
dependence of the spectrum disappears in the large- N limit, at least for sufficiently 
small values of 9 around 9 = 0. In the case of the spectrum, the general large- N 
scaling arguments of Sec. |2], which indicate 9 = 9/N as the relevant Lagrangian 
parameter in the large-iV limit, imply that 0(9 2 ) coefficients should decrease as 
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1/N 2 . The results of Ref. || appear substantially consistent: In the case of the 
string tension they suggest S2 ~ —0.9/N 2 . 

Of course, further investigation is required to put this scenario on a firmer 
ground, using for example other definitions of topological charge. 



6. The case of the two-dimensional CP model 

Issues concerning the 9 dependence can also be discussed in two-dimensional CP^ -1 
models |6B|, |BT[ , 



N - 

2? 



C^—D^zD^z (6.1) 



where z is a iV-component complex scalar field subject to the constraint zz — 1, 
= izd^z is a composite gauge field, and _D M = + iA^ is a covariant deriva- 
tive. They provide an interesting theoretical laboratory. Indeed they present several 
features that hold in QCD: Asymptotic freedom, gauge invariance, existence of a 
confining potential between non gauge invariant states (that is eventually screened 
by the dynamical constituents), and non-trivial topological structure (instantons, 9 
vacua). Moreover, unlike four-dimensional SU(iV) gauge theories, a systematic 1/N 
expansion can be performed around the large- N saddle-point solution |)(| ^7], pSfl . 



Analogously to four-dimensional SU(A r ) gauge theories, one may add a 9 term 
to the Lagrangian, writing 

N 1 

£e = 7T D v z D » z + i9 7T e ^ d » A »' ( 6 - 2 ) 

Zg Z7T 

where q(x) = ^- e^ u d^A v is the the topological charge density. Then one may study 
the 9 dependence of the ground state and other observables. In the following we 
discuss this issue within the 1/N expansion, performed keeping g fixed. Simple 
large- N scaling arguments applied to the Lagrangian (|6.1|) indicate that the relevant 
9 parameter in the large- N limit should be 9 = 9/N . 

Analogously to SU(A^) gauge theories, the ground state energy F{9) depends on 
9. One may define a scaling ground state energy f{9) and expand it around 9 = 0, 

f(0) = M- 2 [F(9) - F(0)] = l -C9 2 [1 + J2 b ^0 2n ) (6-3) 



n=l 



where F(9) is defined as in Eq. (|1.2|) , M is the mass scale at 9 = defined from 
the second moment of the two-point function of the operator Pij{x) = Zi(x)zj(x), 
C is the scaling ratio x/M 2 at 9 = 0, where x is the topological susceptibility. 
The correlation function of the topological charge density, and in particular the 
topological susceptibility, has been computed within the 1/N expansion |7D|, [7T]. 
We have 

C = X /M 2 = — + 0(1/N 2 ) (6.4) 



- 8 - 



The coefficients b 2n are obtained from appropriate 2n-point correlation functions of 
the topological charge density operators at 9 — 0. The analysis of the 1/iV-expansion 
Feynman diagrams |)8| of the connected correlations necessary to compute b 2n shows 
that they are suppressed in the large- N limit, as 

b 2n = 0{l/N 2n ). (6.5) 

This implies that the ground-state energy can be rewritten as 

f(9) = Nf(9 = 9/N), (6.6) 

1- 
-( 

2 



f(e) = \ce 2 (i + J2 h ^ 2n ), 

n=l 



where C = NC and b 2n = N 2n b 2n are 0(N°). Note the analogy with the expected 
9 dependence of the ground-state energy in SU(iV) gauge theories, cf. Eq. (|2.4j ). 
Rather cumbersome calculations lead to the results b 2 = — y, and 6 4 = — =yS. 

Within the 1/N expansion one may also study the dependence of the mass M 
on the parameter 9. We write 

M{9) = M (1 +m 2 9 2 + ...) (6.7) 

The analysis of its diagrams in the corresponding 1/N expansion indicates that m 2 
is suppressed as 

m 2 = 0(1/N 2 ) (6.8) 
Once again, the relevant parameter is seen to be 9 = 9/N . 

7. Critical slowing down of topological modes 

Monte Carlo simulations of critical phenomena in statistical mechanics and of quan- 
tum field theories, such as QCD, in the continuum limit are hampered by the problem 
of critical slowing down (CSD) |72| . The autocorrelation time r, which is related to 
the number of iterations needed to generate a new independent configuration, grows 
with increasing length scale £. In simulations of lattice QCD where the upgrading 
methods are essentially local, it has been observed, see e.g. Refs. [|73|, [74], |75], J|, [7J| 
that the topological modes show autocorrelation times that are typically much larger 
than those of other observables not related to topology, such as Wilson loops and 
their correlators. Actually, the heating method [|4[], used to estimate the topological 
susceptibility, essentially relies on this phenomenon. 

Recent Monte Carlo simulations [|, 7jJ of the four- dimensional SU(iV) lattice 



gauge theories (for N = 3,4, 6) provided evidence of a severe CSD for the topolog- 
ical modes, using a rather standard local overrelaxed upgrading algorithm. Indeed, 
the autocorrelation time r top of the topological charge grows very rapidly with the 
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length scale £ = cr^ 1 ^ 2 , where a is the string tension, showing an apparent expo- 
nential behavior r t0 p ~ exp(c£) in the range of values of £ where data are available. 
Such a phenomenon worsens with increasing N, indeed the constant c appears to 
increase as c oc N. Of course, this behaviour does not depend on the particular 
estimator of the topological charge. This peculiar effect has not been observed in 
plaquette-plaquette or Polyakov line correlations, suggesting an approximate decou- 
pling between topological modes and nontopological ones, such as those determining 
the confining properties. 

These results suggest that the dynamics of the topological modes in Monte Carlo 
simulations is rather different from that of quasi-Gaussian modes. CSD of quasi- 
Gaussian modes for traditional local algorithms, such as standard Metropolis or 
heat bath, is related to an approximate random-walk spread of information around 
the lattice. Thus, the corresponding autocorrelation time r is expected to behave 
as r ~ £ 2 (an independent configuration is obtained when the information travels a 
distance of the order of the correlation length £, and the information is transmitted 
from a given site/link to the nearest neighbors). This guess is correct for Gaussian 
(free field) models; in general one expects that r ~ ^ z where z is a dynamical critical 
exponent, and z m 2 for quasi-Gaussian modes. On the other hand, in the presence of 
relevant topological modes, he random-walk picture may fail, and therefore we may 
have qualitatively different types of CSD. These modes may give rise to sizeable free- 
energy barriers separating different regions of the configuration space. The evolution 
in the configuration space may then present a long-time relaxation due to transitions 
between different topological charge sectors, and the corresponding autocorrelation 
time should behave as r top ~ exp Fb where Fb is the typical free-energy barrier among 
different topological sectors. However, this picture remains rather qualitative, be- 
cause it does not tell us how the typical free-energy barriers scale with the correlation 
length. For example, we may still have a power-law behavior if Fb ~ ln£, or an ex- 
ponential behavior if Fb ~ It is worth mentioning that in physical systems, such 



as random-field Ising systems |77] and glass models [78], the presence of significant 
free-energy barriers in the configuration space causes a very slow dynamics, and an 
effective separation of short-time relaxation within the free-energy basins from long- 
time relaxation related to the transitions between basins. In the case of random-field 
Ising systems the free-energy barrier picture supplemented with scaling arguments 
leads to the prediction that r ~ exp(c£ e ) where 9 is a universal critical exponent 

Q 

The severe CSD experienced by the topological modes under local updating 
algorithms should be a general feature of Monte Carlo simulations of lattice models 
with nontrivial topological properties, since the mechanism behind this phenomenon 
should be similar. This has been also observed in two-dimensional CP^ -1 models 
[ [T9| , p0| . The numerical study of Ref. [[7j| for various values of iV show that an 
exponential Ansatz, i.e. r top ~ exp(c£ 6> ) with 9 « 1/2, and c oc N, provides a 
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good effective description in the range of the correlation length £ where data are 
available (however, the statistical analysis of the data did not allow one to exclude 
an asymptotic power-law behavior r ~ £ z with z ^ N/2 setting in at relatively large 

o. 

The issue of CSD of topological modes is particularly important for lattice QCD, 
because it may pose a serious limitation for numerical studies of physical issues 
related to topological properties, such as the mass and the matrix elements of the rj' 
meson, and in general the physics related to the broken U(1)a symmetry. Indeed, 
it may substantially worsen the cost estimates of the dynamical fermion simulations 



for lattice QCD, see, e.g., Ref. [81 



Finally, we note that although the effects of the topological CSD have not been 
directly observed in plaquette-plaquette or Polyakov line correlations, such a CSD 
will eventually affect them. The point is that the results of Ref. ||, summarized in 
Sec. ||, show that the correlators of plaquette operators and topological charge do 
not vanish at finite N, although they are quite small, and therefore there is not a 
complete decoupling between topological and nontopological modes. Therefore the 
strong critical slowing down that is clearly observed in the topological sector will 
eventually affect also the measurements of nontopological quantities, such as those 
related to the string and glueball spectrum. 
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